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Spin-Parity Analysis of the Centrally produced KsKg system at 800 GeV 
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Results are presented of the spin-parity analysis on a sample of centrally produced mesons in the reaction 
pp — > Psiow{KsKs)pfaat with 800 GeV protons on liquid hydrogen. The spin-parity analysis in the mass region 
between threshold and 1.58 GeV/(? shows that the KsKs system is produced mainly in S wave. The /o(1500) 
is clearly observed in this region. Above 1.58 GeV/c? two solution are possible, one with mainly S wave and 
another with mainly D wave. This ambiguity prevents a unique determination of the spin of the /j(1710) meson. 



1. INTRODUCTION 

The first evidence of the Central Production of 
/o(1500) in the reaction 



pp^Ps{KsKs)pf, Ks 



(1) 



is presented here. The /o(1500) was first observed 
in K^p M interactions and beautifully confirmed 
in low energy pp annihilations by the Crystal Bar- 
rel Collaboration [g[ . Its properties are of current 
interest because it is considered a candidate to be 
the lowest lying glueball state [p| . 

One of the advantages of the final state se- 
lected is that only states with quantum numbers 
J^*^ = (ei;en)++ are allowed to decay into K^Kg. 
This not only greatly simplifies the analysis but 
eliminates confusion coming from all the other 
states. The results presented here are based on 
10% of the 5 X 10^ events recorded by FNAL E690 
during Fermilab's 1991 Fixed Target run. 

2. DATA SELECTION 

2.1. The detector 

The data was taken at Fermilab with an 800 
GeV proton beam on a liquid hydrogen {LH2) 



target, and the E690 spectrometer. The spec- 
trometer is composed of two parts: a) the Main 
Spectrometer (MS), and b) the Beam Spectrom- 
eter (BS). The MS has an approximately conical 
geometrical acceptance with an average 700 mrad 
radius, good momentum resolution from about 
0.2 to 15 GeV/c, a Freon 114 threshold Cherenkov 
counter with a pion threshold of 2.6 GeV/c, a 
time of flight system (TOF) with n/p separation 
up to 1.5 GeV/c, and a target veto system. Nei- 
ther the TOF nor the Cherenkov counter were 
used in the work presented here. 

The BS, used to measure the incoming and out- 
going protons, has an approximately conical ge- 
ometrical acceptance with an average radius of 
about 1 mrad, a pt resolution of 6 MeV/c, and a 
longitudinal momentum resolution of 425 MeV/c. 
The longitudinal momentum acceptance of the 
BS for the interacted beam ranges from approxi- 
mately 650 to 800 GeV/c. 

The trigger required an equal number on in- 
coming and outgoing tracks in the BS and at least 
one additional track in the MS. 



2.2. Event selection 

Final state (|l|) was selected by requiring a pri- 
mary vertex in the LH2 target with two Ks, an 
incoming beam track, and a fast forward proton. 
No direct measurement was made of the slow pro- 
ton ps- The target veto system was used to reject 
events with more than a missing proton. The 
events were accepted when either no veto counter 
was on, or only one veto counter was on with 
the missing pt pointing to it. About 12% of the 
selected events were rejected with the veto sys- 
tem. The missing mass squared seen in Figure n^. a 
shows a clear proton peak with little background; 
the arrows indicate the cuts used in the event se- 
lection. 

The MS has essentially no acceptance for 
XF > 0, which insures a gap of at least 3.5 units 
of rapidity between pf and the central products. 
The average rapidity gap between the KgKs sys- 
tem and Ps is 2.5 units. Figure |l|.b shows the un- 
corrected xp distribution for the KgKg system, 
the arrows indicate the cuts used in the event se- 
lection. 

The TT+TT^ invariant mass for the Ks''s has a 
width oi a — 2 MeV/c^. No direct particle ID 
(Cherenkov or TOF) was used to identify the Kg 
decay products. In about 7% of the events a Kg 
is compatible with a A; these events were kept, to 
avoid any biases in the angular distributions. 

The proton mass was assigned to the missing 
particle in the events that passed the cuts, then 
the three momenta of Ps and the longitudinal mo- 
mentum of Pf were calculated using energy and 
momentum conservation. 

Figure H.c shows the KgKg invariant mass for 
the 11182 events that passed the previous cuts 
and were used in the analysis. The analysis was 
not continued beyond 2 GeV/c^ because the num- 
ber of events is very low. But for —0.22 < xp < 
—0.02 the KgKg invariant mass beyond 2 GeV/c^ 
is smooth, with no evidence of the ^(2230) state 
seen by the BBS Collaboration Q . 

3. DATA ANALYSIS 

3.1. Angular distributions 

The reaction studied here was analyzed as a 
two step process: the production step in which 
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Figure 1. a) Missing mass squared for 

1.4 < M{KsKg) < 1.8 GeV/c^. b) Uncorrected 
Xp distribution, c) Measured KgKg invariant 
mass, d) Acceptance corrected and background 
subtracted KgKg invariant mass. 



an (X) system is formed by the collision of two 
objects (from now on referred to as "pomerons") 
emitted by each of the scattered protons, and the 
decay step in which the object {X) decays into 
KgKg . The production coordinate system was de- 
fined in the CM of the {X) system, with the y-axis 
perpendicular to the plane of the two "pomerons" 
in the overall CM, and the z-axis in the direction 
of the "beam pomeron" in the (X) CM. The two 
variables needed to specify the decay process were 
taken as the polar and azimuthal angles {9, 0) of 
one of the Kg (taken at random) in the produc- 
tion coordinate system. The acceptance corrected 
cos 9 and distributions are shown in Figures y 
and |3|. The acceptance is flat in (j), and dips 
near cos 9 = ±1. On average the correction at 
cos 9 — ±1, relative to the correction at 0, is 65%. 
The five variables used to specify the produc- 
tion process were the transverse momenta of the 




-1 1 

CosO, 1.36<M„<1.42 



75 

50 r 

25 

^ 
-1 











00 
75 




P3 -i°0?i D,a231E-0' 


-, |l ' 


SO 


f^/V^j 


25 


r 





- 1 1 1 1 1 1 1 1 




-1 1 




CosG, 1.6<Mkk<1,66 








uu 


r 


P. oV".^ o^^^Ib 


75 


L 


1 1 1 




Ul^4^^l^-U 


50 


•rViKt 


OR 




^ 


-,,,,!,,,, 



40 
20 - 



n\\ 







5 

*, 1.36<M„<1.42 



60 
40 
20 



« 1 '|l,l,.l""l 



CosO, 1.42<M,K<1,48 



CosQ, 1,66<Mkk<1,72 



5 

*, 1.42<MKK<1,4a 



00 


r 


75 


^tttt)t,tlJ|*t 


50 


\ ''Vtt* / 


25 


--_ \ 



5 

*, 1.6<Mkk<1.66 




1,66<Mkk<1.72 





r 1 


//..rS1.!S / 17 


150 




100 


rt^ 


44iiH/r' 


50 




YT)" ' 



-1 1 

CosO, 1.48<Mkk<1.54 



' *=/ndf 



100 - 




-1 1 

CosG, 1.54<M>,<1.6 



100 

75 

50 

25 





-1 1 

CosG, 1.72<M„<1.78 




-1 1 

CosO, 1.7a<M«,<1.84 



Figure 2. Acceptance corrected cos 9 angular dis- 
tributions in bins of the K^Kg invariant mass, 
starting at 1.36 GeV/c^ in steps of 60 MeV/c^ 
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Figure 3. Acceptance corrected cj) angular dis- 
tributions in bins of the K^Kg invariant mass, 
starting at 1.36 GeV/c^ in steps of 60 MeV/c^. 



slow and fast protons {pf g,pf r), the xp and in- KgKg CM. Although our 11182 events constitute 

variant mass of the K^Kg system, and 5, the angle a large sample, it is not large enough to bin the 

between the planes of the scattered protons in the data in all five production variables. The present 



analysis was done in bins of the K^Ks invariant 
mass for the xp selected region, and integrating 
over p1 g, p1 r and 5. 

3.2. Partial wave analysis 

The acceptance corrected moments, defined as 
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are shown in Figure 0. The odd moments (not 
shown) are consistent with zero, as expected for 
a system of two identical bosons. The accep- 
tance corrected mass distribution (ioo moment) 
is shown in Figure 0.d. The error bars are statis- 
tical errors only. 

In the two step process considered here the 
{X) system is formed by the interchange of two 
"pomerons" and it decays afterwards indepen- 
dently of the two final state protons. The two 
"pomerons" form a plane; parity in the strong 
interactions implies that reflection in this plane 
should be a symmetry of the system [g . There- 
fore the amplitudes used for the Partial Wave 
Analysis (PWA) were defined in the reflectivity 
basis [g| . Since the tis and ^44 moments are con- 
sistent with zero (see Fig ^) , only spherical har- 
monics with I = 0,2 and m = 0, ±1 were consid- 
ered. The waves used were L^„, with L — S,D, 
m > and e = ±1: 
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Figure 4. Acceptance corrected moments as a 



Waves with different reflectivity e do not interfere. 



The PWA analysis was done in two different 
ways. First since the angular distributions are 
fairly flat only Sq and Dq waves were used: a) 
by fitting to the cos 9 angular distributions, and 
b) by using the extended maximum likelihood 



method. The results of the fit to the cos 9 angular 
distributions are shown in Figure 0. Within errors 
the results were the same in both cases, giving a 
solution that, except for two small D wave contri- 
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Figure 5. Waves as a function of KgKg invariant 
mass for solution one. a) S and b) total D waves, 
c) to e) individual D wave, and f) and g) phases 
relative to the S wave. 
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Figure 6. Waves as a function of KgKs invariant 
mass for solution two. a) 5* and b) total D waves, 
c) to e) individual D wave, and f) and g) phases 
relative to the S wave. 



butions at c^l.3 GeV/c^ and ~1.6 GeV/c^, was all 
S wave. Second, all four waves (&|[) were used. 



The amplitudes were extracted both a) from the 
moments shown in Figure 0, and b) by maxi- 



mizing the extended likelihood with respect to 
the four wave moduli and the two relative phases 
(P{Dq^) — (p{Sq). Within errors both analyses 
gave the same answer. 

When using the four waves (|^||) the inherent 
ambiguities of a two body system are such that 
there are two solutions for each mass bin [gJ7|. 
Both solutions give identical moments or identi- 
cal values of the Likelihood. In order to continue 
the solutions from one mass bin to the next, one 
follows the Barrelet zeros. In general these zeros 
are complex and one lies above the real axis and 
the other lies below it. When the zeros cross the 
real axis the solutions bifurcate ||J^ . In the anal- 
ysis presented here, there is a bifurcation point at 
1.58 GeV/c^. Before this bifurcation point there 
are only two solutions, one which is mostly S 
wave, and another that is mostly D wave. Since 
at threshold the KgKg cross section is dominated 
by the presence of the /o(980) ||] it is possible 
to eliminate the solution that has a very small S 
wave contribution at threshold. The remaining 
solution bifurcates at 1.58 GeV/c^ into a solution 
that has a large S wave contribution (solution 
one), and another that has a large D wave com- 
ponent (solution two). The solutions obtained us- 
ing maximum likelihood are shown in Figures ^ 
and |6|. Solution one is shown in Figure |^, and 
solution two in Figure g. The errors shown are 
statistical errors only. 

A striking feature of both solutions is the large 
S wave peak observed at ~1.5 GeV/c^. This cor- 
responds to the /o(1500) observed by the Crys- 
tal Barrel collaboration M. The mass peaks at 
1.52 GeV/c^ instead of at 1.50 GeV/c^ but this 
could be easily due to interference with the S 
wave background. 

Beyond 1.58 GeV/c^ both solution one and two 
are equally valid, and at the moment there is no 
way to decide with this data alone which of the 
two solutions is the correct one. However, given 
that beyond 1.58 GeV/c^ the angular distribu- 
tions are fairly structureless, and that an analy- 
sis in cos^ alone gives very little D wave, solution 
one could be favored. 



4. CONCLUSIONS 

A PWA analysis in a sample of 11182 cen- 
trally produced KgKs events at 800 GeV has 
been presented. Two solutions have been found 
in the analysis. In both of them a clear /o(1500) 
has been observed. The ambiguity above 1.58 
GeV/c^ prevents a unique determination of the 
spin of the /j(1710) meson. Due to lack of statis- 
tics the analysis was not carried out beyond 2 
GeV/c^, but the KgKs invariant mass spectrum 
is smooth beyond that point and shows no sign 
of the ^(2230) meson seen by the BES Collabora- 
tion [|. 
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